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ABSTRACT 

Two-phase gas-liquid flows occur in a wide variety 
of situations. In addition to normal gravity 
applications, such flows may occur in space operations 
such as active thermal control systems, power cycles, 
and storage and transfer of cryogenic fluids. 

Various flow patterns exhibiting characteristic 
spatial and temporal distribution of the two phases are 
observed in two-phase flows. The magnitude and 
orientation of gravity with respect to the flow has a 
strong impact on the flow patterns observed and on 
their boundaries. The identification of the flow pattern 
of a flow is somewhat subjective. The same two-phase 
flow (especially near a flow pattern transition 
boundary) may be categorized differently by different 
researchers. 

Two-phase flow patterns are somewhat simplified 
in microgravity, where only three flow patterns 
(bubble, slug and annular) have been observed. 
Annular flow is obtained for a wide range of gas and 
liquid flow rates, and it is expected to occur in many 
situations under microgravity conditions. Slug flow 
needs to be avoided, because vibrations caused by 
slugs result in unwanted accelerations. Therefore, it is 
important to be able to accurately predict the flow 
pattern which exists under given operating conditions. 

It is known that the wavy liquid Film in annular 
flow has a profound influence on the transfer of 
momentum and heat between the phases. Thus, an 
understanding of the characteristics of the wavy film is 
essential for developing accurate correlations. 

In this work, we review our recent results on flow 
pattern transitions and wavy films in microgravity. 

FLOW PATTERN TRANSITIONS IN 
MICROGRAVITY 

The absence of gravity reduces the number of 
observed flow patterns to three and it also reduces the 
number of dimensionless groups needed to 
characterize two-phase flows by two. For 
microgravity two-phase flows in a smooth pipe, there 
are five relevant dimensionless groups: 
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Here, Uls is the superficial velocity of the liquid based 


on a single phase flow and U G s is the superficial gas 
velocity. The viscosity, density and surface tension 
are denoted by [i, p and a, respectively. When the last 
two ratios are small, the flow pattern transition 
boundary may be assumed to be a weak function of 
these two ratios. The number of dimensionless groups 
may be further reduced to two by combining Rc LS and 

We LS and ignoring the influence of We LS alone. 

For microgravity two-phase flows, (using 
experimental data from various systems) we developed 
a pair of dimensionless flow pattern transition maps 
shown in Figs. 1 and 2 (Jayawardena et al, 1997). 
These maps suggest the importance of Suratman 
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transitions between the flow patterns. 


BUBBLE-SLUG TRANSITION 

The bubble-slug flow pattern transition occurs at a 
particular value of the ratio (Recs^ e Ls)* This 
transitional value, (Re GS /Rc LS ) t , depends on the 
Suratman number as given below: 

(Rcgs/Rcls), = K, Su M (1) 

SLUG- ANNULAR TRANSITION 

It is found that when the Suratman number is less 
than 10 6 , slug-annular flow pattern transition occurs at 
a particular value of the ratio (Re G s/R e Ls)’ This 
transitional value, (RcGs^ e Ls)f was f°und to depend 
on the Suratman number as given below: 

(R^GS^ e Ls)t = ^2 2/3 • 

When the Suratman number is greater than 10 6 , slug- 
annular flow pattern transition was found to be a 
function of the gas Reynolds number, Re G s- This 
transitional value, (Re G s)t» was found to depend on the 
Suratman number as given below: 

(Re os ) t = K 3 Su 2 . (3) 

Experimental data shown in Figs. 1 and 2 suggest 
the following numerical values: 
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Ki =464.16, 

(4a) 

K 2 = 4641.6, 

(4b) 

K 3 = 2x1 0‘ 9 . 

(4c) 


The Suratman number is determined by the tube 
diameter and physical properties of the fluid. Thus, the 


proposed maps can be used to identify the flow pattern 
for any given two-phase system, even when there are 
no prior experimental microgravity flow pattern data. 



Figure 1. Dimensionless flow pattern map for microgravity two-phase flows. (For slug-annular transition for Su > 10 6 and for 
bubble-slug transition). 
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Figure 3. Results of recent microgravity two phase flow tests conducted using various fluids as shown in table 1 . 


8 





STUDIES ON NORMAL AND MICROGRAVITY ANNULAR TWO PHASE FLOWS: 
V. Balakotaiah, S. S. Jayawardena and L. T. Nguyen 


VALIDATION OF PROPOSED MAP 

We have collected new flow pattern data on 
micro-gravity two phase flows. These tests were done 
using the same test rig used by Bousman (1994), but 
with different working fluids. The Suratman number 
was changed by reducing the surface tension and/or 
by increasing the viscosity using a surfactant and/or 
glycerin. Figure 3 shows these new data points. The 
fluid properties and corresponding Suratman numbers 
are listed in Table 1 . 

These new experiments confirmed that the slug- 
annular transition boundary given by Eq. (2) extends 
to low values of the Suratman number (10 2 < Su < 
10 4 ). However, they did not verify the bubble-slug 
transition boundary. The proposed boundary may not 
applicable for low Suratman number bubble-slug 
transitions. On the other hand, due to experimental 
limitations, the flow visualization was done after 85 
pipe diameters. Since the feed section generated 
annular flow. It is possible that the flow development 
length was not sufficient to distribute the phases to 
achieve bubble flow. 

When we planned these experiments, our aim was 
to validate the slug-annular transition for various two 
phase flows in microgravity and we could do that with 
flows having Suratman numbers as low as 770. 


Table 1 . Properties of the fluids used in these tests. 


Test Fluid 
(Wt % of 
Glycerin) 

Su 

Properties 1 

Pl 

Pl 

G 

AWG 

(70) 

1800 

1181 

23.00 

64.6 

AWGZ 

(70) 

770 

1181 

23.00 

27.0 

AWG 

(60) 

8900 

1153 

10.66 

68.9 

AWGZ 

(50) 

11000 

1125 

6.03 

27.0 

AWG 

(24) 

250000 

1056 

1.98 

72.2 

AWG 

(64) 

5400 

1164 

13.63 

68.4 

AWG 

(32) 

140000 

1077 

2.63 

71.8 

AWGZ 

(32) 

56000 

1077 

2.63 

28.5 


(weight percentage of glycerin is given in parentheses) 
Key: AW = air- water, 


AWG = air-water-glycerin, 

AWZ = air-watcr-Zonyl, 

AWGZ = air-waler-glycerin-Zonyl, 
Units: p L (kg/m 3 ), p L (cP) and a (dyne/cm), 


SYSTEM SPECIFIC. DIMENSIONAL FLOW 
PATTERN MAPS 

Some system specific flow pattern maps (for a 
given fluid in a selected tube size) are plotted using the 
two-phase quality vs. mass flow rate. Such maps are 
useful for single-component, two-phase system. Our 
flow pattern boundaries may be used to create such 
system specific flow pattern maps. 

The two-phase quality x, is defined as 

(Pg^gs) 

(Pg^gs + Pl^ls) 

Combining Eqs. (1), (5) and material balance on the 
two-phase flow gives the following expressions for the 
bubble-slug flow pattern boundary for any Suratman 
number: 
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This transitional quality, x tBS , is independent of the 
total mass flow rate. 

The criteria for the slug-annular transition depends 
on the Suratman number of the system. The following 
expression (from Eqs. (2,) (5) and material balance) 
gives the slug-annular flow pattern boundary when the 
Suratman number is smaller than 10 6 : 
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- (Su<10 6 ) (7)) 

) 

, is also independent of 


the total mass flow rale. 

Combining Eqs. (3), (5) and material balance gives 
the following slug-annular flow pattern boundary 

when the Suratman number is larger than 10 6 : 

3 2 2 

® x t,s-A= K 3 i T"^F P (Su > 10 6 ) (8) 


where m is the total mass flow rate of the two- 
phase flow. Equation (8) suggests that it is possible to 
obtain slug flow at higher values of quality, provided 
mass flow rate of the two-phase flow is low. 

An example of such a single-component, two-phase 
system is a refrigerant flowing in a tube. In such 
cases, the total mass flow rate and the quality are the 
preferred variables to prepare flow pattern maps. 
Figure 4 shows the flow pattern boundaries on a two- 
phase quality vs. mass flow rate map, for a selected 
two-phase system in microgravity. 

The bubble flow pattern is observed in an 
extremely small range of quality. The value for x t B s 
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at that transition is about 0.04%. For most flow 
conditions in microgravity, annular flow is observed. 



Figure 4. Dimensional, system specific flow pattern map for 
a single component microgravity two-phase flow. (For 
refrigerant R134a at 70°F in 25.4 mm ID tube, Su = 
4.5x1 0 6 ) 

The system specific flow pattern map for a low 
(<10 6 ) Suralman number system consists of two 
straight lines, since the two-phase quality at the flow 
pattern boundary is independent of the total mass flow 
rate. A detailed discussion including the critical void 
fraction of a two component systems is given in 
Jayawardena and Balakotaiah (1997). 

WAVY FILMS IN MICROGRAVITY 

We now study characteristics of the wavy films in 
microgravity annular flows. Bousman (1994) collected 
data on wavy film profiles in microgravity using air- 
water, air-water-glycerin and air-water-Zonyl systems. 
His film thickness data are shown in Fig. 6. Where the 
dimensionless film thickness is defined as 



Since wall measured shear stress (t w ) is not accurate, 
we estimate T w using measured pressure gradient: 



where x w is the average interfacial stress and D is the 
tube diameter. The liquid Reynolds number Ret is 
defined as 



where T is the flow rate in the film per unit perimeter. 
Since film flow rate is not measured, we approximate 
Re L by Re^, assuming negligible entrainment. 

A relationship between h + and Re L is derived by 
assuming a velocity profile. The following velocity 
profile is assumed in the film for laminar flow 



For the turbulent flow, the universal velocity profile 
u + = y + 0 < y* < 1 0 

u + = 2.78 In y + + 3.60 10<y + (13) 

where, u + is the dimensionless friction velocity in the 

film (u + = u/u*) is assumed in the liquid film. 

These equations gives 

Re t = 2 ( h + ) 2 0 < h + < 10 

Re. = 3.28 h + + 11.12 In h + - 88.8 10<h 

(14) 

A plot of this equation is shown in Fig. 5 as the 
solid line. The agreement with data is good. The mean 
film thickness of normal and microgravity two phase 
flows are not fundamentally different. 

Since we expect the surface tension effects to 
dominate under microgravity conditions, it is of 
interest to compare the rms values of the film thickness 
under normal and microgravity conditions. (Note: 
The rms film thickness is actually the standard 
deviation since we subtract the mean film thickness. 
However, we are following here the literature 
terminology.) 

The rms values of the film thickness fluctuations is 
normalized using the same length scale defined using 

the friction velocity u*, and it is denoted as h*^. 

When h^s is plotted against Re^s. it shows a large 
amount of scatter, partly due to the fact that these data 
points were collected at different gas velocities ranging 
from about 5 m/s to about 25 m/s. Experimental 
observations show that the annular flow liquid film 
becomes smoother at higher gas flow rates. 

The following correlation predicts measured h + rms 
values for various flow conditions used in these tests: 
h + ms = ARe LS 07 (15) 

Here A is a coefficient which depends only on the 
gas flow rate. The available data suggests: 

A = (0.135-0.0033 U GS ). (16) 

Figure 6 is a comparison of the predicted h + rms 
values with the experimental values. One dilemma 
here is that our intuition suggests that h + rms should be 
affected by the surface tension, where as Eq. (15) 
indicates the contrary. 
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Another quantity of interest is the enhancement in 
the friction factor. Here we compare the friction factor 
based on the interfacial stress with that of a 
corresponding single phase gas flow. The friction 
factor is defined as 

<n > 

2 Pg U gs 

Asali et al (1985) have analyzed data for co- 
current down flow and presented a correlation of the 
form 



Figure 5. Dependence of normalized liquid film thickness on 
the liquid Reynolds number for microgravity annular two- 
phase flow and the comparison with the proposed 
correlations. 

j - 1 ) Re 0 2 GS = 0.45 (h + G - 4).. (18) 

Here, f s is the gas phase friction factor that would 
exist in a smooth tube and 



The microgravity data is shown in Fig. 7 along with 
the correlation given by Eq. (20). Though the data is 
more scattered, it is clear that the enhancement in the 
interfacial friction factor is higher in microgravity than 
in normal gravity. 

If we use all the available micro-gravity data, then 
the interfacial shear stress can be represented by the 
following equation: 

f Re°'^ GS = 0.42 h + G 1 18 (20) 


When we eliminated data points close to the slug 
annular transition, the following two correlations are 
obtained. 

^-ljRe°- 2 GS = 0.16h + G 13 (21) 

for water and water-glycerin mixtures and 

^-ljRe°- 2 GS = 0.25h + G 13 (22) 

for waler-Zonyl mixture. (Zonyl reduces the 
surface tension by about a factor three). That indicates 
that surface tension influences the interfacial structure 
of the waves. 


o (AW) a (AWG) * (AWZ) 
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Figure 6. Comparison of calculated h ^ values with the 
experimental data. 
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Figure 7. Comparison of the friction factor for the interfacial 
shear stress in microgravity annular flows with the existing 
normal gravity annular flow correlation of Asali et al. (1985) 
(Dashed line) with the new correlation (solid line). 

CONCLUSIONS AND DISCUSSION 

The recently proposed flow pattern map was 
verified at low Suratman number systems. 
Experiments confirmed the applicability of the 
proposed slug-annular flow pattern transition boundary 
for two phase flow systems with a low Suratman 
number. 

System specific, dimensional flow pattern 
boundaries can be obtained for microgravity two phase 
Hows. For a wide range of two-phase qualities, slug 
flow can be observed in such a system if the total mass 
flow rate is small, provided Suratman number is large. 
Bubbly flow exists in microgravily only at extremely 
low two-phase qualities. 

The dimensionless parameter affecting the slug- 
annular transition changes at a Suratman number of 
about 10 6 . Thus, we expect a dramatic change in the 
mechanism governing the slug-annular transition. It is 
very important to study liquid film characteristics, 
pressure gradients and drop entrainment measurements 
for high Suratman number two-phase flows in 
microgravity. 

When properly non-dimen sionalizxd using the 
pressure gradient, the mean film thickness of an 
annular flow depends only on the liquid Reynolds 
number, irrespective of gravity. Agreement between 
experimental data and the prediction suggests that 
entrainment of liquid droplets in microgravity annular 
flows to be small. This speculation needs 
experimental validation. 


The differences between the normal and 
microgravity annular liquid films arc apparent when 
we compare the interfacial friction factors. It is found 
that the friction factor enhancement is higher in 
microgravity conditions than in normal gravity 
downward annular flows. In microgravity annular 
flows, the enhancement to the friction factor is surface 
tension dependent. 
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